Mon. Not. R. Astron. Soc. 0Q0.[TinT1(2009) Printed 17 September 2009 (MN KTeX style file v2.2) 



A new cosmic shear function: Optimised E-/B-mode decomposition 
on a finite interval 



Liping Fu 1 ' 2 * and Martin Kilbinger 3 



O 

o 

Oh 

in 



O 

U 

6 

_fH_ 

(N 

> 
in 

On 

o 
o 

On 
O 



13 



INAF, Osservatorio Astronomico di Capodimonte, via Moiariello 16, 80131 Napoli, Italy 
2 Shanghai Key Lab for Astrophysics, Shanghai Normal University, Shanghai 200234, P. R. China 
3 Institut d Astrophysique de Paris, 98bis boulevard Arago, F-75014 Paris, France 



Accepted / Received 



ABSTRACT 

The decomposition of the cosmic shear field into E- and B-mode is an important diag- 
nostic in weak gravitational lensing. However, commonly used techniques to perform this 
separation suffer from mode-mixing on very small or very large scales. We introduce a new 
E-/B-mode decomposition of the cosmic shear two-point correlation on a finite interval. This 
new statistic is optimised for cosmological applications, by maximising the signal-to-noise 
ratio (S/N) and a figure of merit (FoM) based on the Fisher matrix of the cosmological param- 
eters £2 m and cr 8 . 

We improve both S/N and FoM results substantially with respect to the recently intro- 
duced ring statistic, which also provides E-/B-mode separation on a finite angular range. The 
S/N (FoM) is larger by a factor of three (two) on angular scales between 1 and 220 arc min- 
utes. In addition, it yields better results than for the aperture-mass dispersion (MjL), with 
improvements of 20% (10%) for S/N (FoM). Our results depend on the survey parameters, 
most importantly on the covariance of the two-point shear correlation function. Although 
we assume parameters according to the CFHTLS-Wide survey, our method and optimisation 
scheme can be applied easily to any given survey settings and observing parameters. Arbitrary 
quantities, with respect to which the E-/B-mode filter is optimised, can be defined, therefore 
generalising the aim and context of the new shear statistic. 
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1 INTRODUCTION 

Cosmic shear, the weak gravitational lensing effect induced on im- 
ages of distant galaxies by the large-scale structure in the Uni- 
vers e, has become a standard tool for observational cosmol ogy 
(see lSchneiderl2006l : lHoekstra & Jain 2008: lMunshi et al]2008l for 
recent reviews). Large surveys have used cosmic shear to obtain 
measurements of the matter density £l m and the density fluctua- 
tion amplitude <x g . Recent constrai nts were obtained from ground- 
based surveys such as CFHTLSpl dBeniamin et ai1l2007l ; I Fu et al.1 
120081) and GaBoDS^l ( iHetterscheidt et alj|2007t). Space-based sur - 
veys like COSMO^I jLeauthaud et alj|2007l; iMassev et alj|2007l) 
and parallel ACS data ( Schrabback et al. 2007b advanced cosmic 
shear observations to very small angular sca les. Cosmic shear 
has contributed to c onstraining dark energy jjarvis et alj 120061 : 
Kilbin ger et al.l200^) . It is considered to be one of the most promis- 
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ing method to shed light onto the origin of the recent accelerated ex - 
pansion of the Universe JAlbrecht et a l. 2006; Pe acock et al.1 2006), 
and is a major science driver for many future surveys like KIDS, 
Pan-STARRS, DES, LSST, IDEM or Euclid. 

One of the (few) diagnostics for cosmic shear analyses is the 
decomposition of the shear field into its E- and B-mode. Most com- 
monly, the shear power spectrum or, equivalently, the shear correla- 
tion function, is split into the gradient (E) an d curl (B) component 
dCrittenden et alj|200_t ISchneider et alj|2002h . Gravitational lens- 
ing produces, to first order, a curl-free shear field and therefore, 
the presence of a B-mode is an indication of residual systematics 
in the PSF correction and shape measurement analysis. To obtain 
competitive constraints on cosmological parameters, in particular 
for dark-energy or beyond-standard physics, galaxy shapes have to 
be determined to sub-percentage precision. This requires excellent 
correction of PSF effects arising from the atmosphere, telescope 
and camera imperfections. 

Apart from observational effects and measurement systemat- 
ics, the cosmic shear signal can be severely contaminated by in- 
trinsic correlations of galaxy orientati on with other galax ies or 
their surrounding dark matter structures faeavens et alj200Ch . This 
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occurs for galaxies at the same redshift e.g. which reside in the 
same dark halo (intrinsic alignment). It also affects galaxies at very 
different redshifts, where a background galaxy is lensed by mat- 
ter surrounding a foreground galaxy, therefore inducing a shape- 
shear-correlation between the two galaxies . An E-mode as well as 
a B-mode arises fr om intrinsic alignment dCrittenden et alj|2002l ; 
iMackev et ail2002l) . 

Standard methods to separate the E- and B-mode power spec- 
tra (or correlation functio ns) involve integrals up t o arbitrary small 
or la r ge angular scales I Schneider et alj 1 19981 ; ICrittenden et al.l 
120021 ; [Schneider et alj|2002h . However, shear correlations can be 
observed only on a finite interval. Since the shear field is probed 
at galaxy positions, the smallest usable scale is given by the confu- 
sion limit for close galaxy pairs, which is typically several arc min- 
utes for ground-based surveys. The largest observed distance is for 
current surveys several degrees. These limits make the E-/B-mode 
separation imperfect, and a mix ing of modes is induced at the 1- 
10% level dKilbinger et al.112006) . To circumvent this shortcoming, 
a new second-order function, the so-called "ring statistic", was in- 
troduc ed which permits a clear E-/B -mode separation on a finite in- 
terval dSchneider & Kilbingeij|2007l, hereafter SK07). In addition, 
the authors developed conditions for general filter functions neces- 
sary for an E-/B-mode decomposition for a finite angular range. 

In this work, we present a method to find filter functions which 
fulfill the SK07 conditions. We devise a scheme which provides an 
optimised E-/B-mode decomposition on a finite interval. The opti- 
misation is performed with regard to cosmological applications of 
cosmic shear; the signal-to-noise ratio and a Fisher matrix figure or 
merit are the quantities to be maximised. This paper is organised as 
follows: In Sect. [2] we briefly review the results from SK07 before 
we present our optimisation method. Results for the signal-to-noise 
ratio and the figure of merit are shown in Sect. [3] We conclude the 
paper with a summary (Sect.[4]l and an outlook (Sect. [5). 



2 METHOD 

2.1 E-and B-mode Decomposition of the shear correlation 
function on a finite interval 

We define the general second-order cosmic shear functions K E and 
Kb. 

R b = J f d*0[r + (0)f + (*) + r_(0)f_(#)]; 

* Jo 

* Jo 



(1) 



as integ rals over the s hear two-point correlation functions and 
(e.g. iKaiserl 1 1992j) with arbitrary filter functions T + and T_. 
These expressions correspond to Eq. (39) from SK07, with an ad- 
ditional factor 1/2 in our definition. In terms of the E- and B- 
mode power spectrum, P E and P B , respectively, the shear two-point 
correlation function ar e given as the following Hankel transforms 
dSchneider et alj2002h 
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with J v being the v th -order Bessel function of the first kind. Insert- 
ing Eq. (O into Eq. QJ yields 
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and an analogous expression for 7?b- The Hankel transforms of T + 
and T_ are defined as 



di?#r + ,-0?)Jo.4(i?£). 
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To provide an E- and B-mode decomposition, in the sense 
that K E only depends on the E-mode of the shear field and 'Kb 
only on its B-mode, the two Hankel transform have to be identi- 
cal, r W E = "W B . After some algebra, one finds that the following 
equiv alent relations betwee n the filter functions T + and T. must 
hold dSchneider et~aT]|2002h 
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Therefore, for an arbitrary function T + , a corresponding filter T_ 
can be derived from T + to provide an E- and B-mode decomposi- 
tion, and vice versa. In the absence of a B-mode we have K B = 0, 
and *R E can be obtained from £ + or alone, 
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We further require K E and K B to depend on the shear correla- 
tion given at angular scales ■& in a finite interval, < ffmin < & < 
j? max < 00. Thus, we demand T_ to have finite support [i? m i n ; )? max ] 
and T + (i?) to vanish for 1? < # m i„, Eq. lf5J then implies the following 
integral constraints on the filter function T_ (SK07): 
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r_(i?) = 0. 
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In addition, it follows that T + (i?) = for ■& > # max . Using the finite 
support of T + in (Eq.[6}, we get integral constraints for T+, 



J#mi 



d$$T + (&) 



di?i? 3 r + (!?) = 0. 



(9) 



Then, *R E and 7? B are functions of the two angular scales # min 
and !? max . We will discuss the scale-dependence in more detail in 
SectM6l 

SK07 constructed a set of functions, Z + ,Z_ in their nota- 
tion, which satisfy Eqs. l[8] |9). Those functions were motivated 
from a geometrical ansatz, by considering two concentric, non- 
overlapping rings. If the shear correlation is calculated from galaxy 
pairs of which one galaxy lies in the inner ring and the other galaxy 
in the outer ring, the E-/B-mode decomposition on a finite interval 
is guaranteed by construction. The form of the function Z + origi- 
nated in a specific choice of the weight profile over the two rings. 
The relation between T ± and Z ± is 



r ± (#) = ir 2 z ± o?/i? max ). 



(10) 



Note that in SK07 the analogous integrals to Eq. (fTJ using Z+ are 
carried out over the integration variable i?/!? max and extend from 
& m m/&mix t° !• The shear second-order functions corresponding to 
Eq. (0 are denoted as (7?'??)e and (7?'/?)b, respectively. 



There are infinitely many functions which fulfill the above in- 
tegral constraints. Their choice can of course be detached from the 
geometrical considerations of the "ring statistic". In this paper we 
define a general filter function which we will optimise regarding 
some specific criterion. This criterion will be related to the cosmo- 
logical information output from a cosmic shear survey. We will use 
two cases, the signal-to-noise ratio and the Fisher matrix of cosmo- 
logical parameters. The results are presented in Sect.[3] 



2.2 Parametrisation of the filter function 

For the optimisation problem, we focus on T + since T- can be de- 
rived from T+ (Eq.[6}. First, we remap T + to the interval [-1; +1] 
by defining 

f+(x) = T + (Ax + B) = T + (&) for x e [-1; 1]; 

A=(# max -# mm )/2; B = (# max ' u mm )/2. (11) 
With that the two integral constraints (Eq.[9]» become 

dx (x + R) f + (x) = J dx (x + Rf f + (x) = 0, (12) 



where we have defined the ratio R as 
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Next, we decompose T+ into a finite sum of orthogonal poly- 
nomials: 



f+(x) = ^ a„ C„(x). 



(14) 



This representation allows us to find an optimal filter function by 
varying the coefficients a„. Sect. 12.71 The polynomials C„ can be 
chosen freely; we use Chebyshev polynomials of the second kind, 



U„(X) : 



sin[(n + 1) arccos x] 
sin(arccos x) 



(15) 



The optimisation process is then performed by varying the coeffi- 
cients a„; this is described in detail in Sect. 12.71 

Apart from the integral constraints (Eq.[9]l, one could require 
T + to be zero at the interval boundaries, 



Continuity : 7V(# raax ) = T+(& mm ) = 0. 



(16) 



Additional constraints could be added, for example differentiability 
at the boundaries. We will discuss their effects on the results in 
Sect.l3Tl 



2.3 Satisfying the constraints 

In general, the function T + , or equivalently T+, is constrained by 
K > 2 equations in the form F,„[f + ] = 0, m = . . .K — 1. If the 
functionals F,„ are linear in T + , applying Eq. l !14t leads to 



^ j fmn @n 0* fmn ■ F m \C }1 \. 
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The matrix element /„„, is the m lh constraint applied to the orthog- 
onal polynomial of order n. For example, taking the first constraint 
in Eq. l !12t we get 



./()/, 



dx (x + R)C n (x), 



(18) 



which can be integrated analytically. The other matrix elements f„,„ 
are obtained analogously. 

K constraints fix K coefficients of the decomposition fEq. U4b . 
the remaining N - K coefficients can be chosen arbitrarily. We use 
the highest N - K coefficients as free parameters (n = K . . . N — 1) 
and fix the first K coefficients (n = 0. . .K — 1) as follows. Define 



Eq. dl7t can then be written as 

K- I 
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m = 0...K-l, 



m = 0...K-l. 



(19) 
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This (K x K)-mairix equation is solved for the first K coefficients 
dn (n = . . . K — 1) by inverting the square (sub-)matrix {f„ m ) m ,n<K 
on the left-hand side of the above equation. If the constraints are 
chosen such that they are linearly independent (which is this case), 
this matrix equation has a unique, non-trivial solution. 

In addition to the K constraints, we impose an integral nor- 
malisation of the filter function given by the L 2 -norm 



iif + ig = J[ 



dxw(x)f^(x) = 1, 



(21) 



where w is the corresponding weight of the polynomial family, 
vv(x) = (1 - x 2 )' /2 in the case of second-kind Chebyshev poly- 
nomials. This normalisation does not affect the constraints which 
are independent of a multiplication of all a„ with a common factor. 
Note also that the quantities which we will optimise in Sect. [3] do 
not depend on the normalisation. 

2.4 Calculation of T- from T + 

To obtain the function ?_ := T_(Ax + B) we transform Eq. ((6) to 



r x , x' + r 
r_(x) = r + (x) + 4 dx'r + (x')- — — 

J_, (x + Rf 
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x = -l..A, 
which can be written as 

N-l 

T-{x) = ^ a n [C„(x) + a„(x)] ; 

n=0 

X x x' + R 
dx'CJx')- - 
1 (x + R) 

inserting the decomposition (Eq.ll4t. We define 

F< v) (x) = J dx' (x')"C„(x'); v = 0,1,2,3, 

and write Eq. d24b as 

a n (x) =y[R(r - R)F®> + (1 - 2rR - rR 2 )F? 

-rR(R + 2)Ff ] - rF< 3) ] ; 
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where we dropped the argument x from F„ v) and r. The integral 
(Eq.|25j for v = and C„ = U„ is 

( (-l)" + xT„(x)-(l-x 2 )U„_ 1 (x) 
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where T„ is n th -order Chebyshev polynomial of the first kind, 
T„(x) = cos(rc arccos x). 



(28) 



By using the recurrence relation of the Chebyshev polynomi- 
als, we obtain the other functions, 



(#min> #max) — > -t> two ways to handle the scale-dependence seem to 
be suitable: (1) Leaving one scale constant, and varying the other 
scale, e.g. & min = const., A = i? max . (2) Leaving the ratio of both 
scales constant, A = rj. We will pursue both ways in this paper. 



F?=\[Ff} 2 + 2^ + F%}- (29) 
F< 3 > =1 [f ,0 > + 3F (0 > + 3F (0) , + F <0 > 1 . 

« o [ n+3 n+1 n— 1 H-3J 

Note that Eq. d27| > is valid for integer n, since the expressions for 
the orthogonal polynomials are well-defined for negative n. 

With that, it can be readily checked whether the coefficients a„ 
obtained with the method described in Sect. 12.31 and the resulting 
filter functions T ± indeed provide an E-/B-mode decomposition, by 
verifying that the B-mode H% (Eq.[TJ is zero. In Sect. |3.6| we discuss 
numerical issues when calculating the B-mode. 

2.5 Relation to the lensing power spectrum 

Eq. (|3) shows the relation of R E to the power spectrum. Assuming 
a pure E-mode, this equation reads 



2.6.1 Fixed minimum scale, $ m 
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where the Fourier-space filter function r W E can be written as 

AM 

max 



(30) 



(3D 



<W IP S . S ) 



!' 



d&#C n [(fi-B)/A]J (flQ, (32) 



with A and B given in Eq. i ll It . Analogously, R% can be written in 
terms of the B-mode power spectrum F B 



R n 



2nJ 



(33) 



with 1V E = W B , see Sect.l2Tl 

Unfortunately there is no simple analytical expression of 
Eq. d32t . which would be d esirable to calcu late Eq. J30b using a 
fast Hankel transform (FHT, Hamilton 2000). For the moment, the 
most efficient method to calculate Re from a model power spectrum 
is to obtain £ + from P E by FHT and to integrate it via Eq. Q which 
is fast since in our case T + is a rather low-order polynomial, as we 
will see later. 



2.6 Angular-scale-dependence 

The two constraints (Eq. O depend on the ratio of angular scales 
f] = ^min/^niajf A filter function f + with given coefficients a„ sat- 
isfying those constraints does not in general fulfill the same con- 
straints for a different n. This therefore causes the inconvenience of 
having a different filter function for each angular scale. 

It should be noted that although formally Re is a function of 
two angular scales, the information about cosmology and large- 
scale structure is captured by a single parameter, let us call it A. 
This is because the large-scale matter power spectrum, of which 
R E is a logarithmic convolution (Sect. [23t , only depends on a sin- 
gle scalar. We therefore expect a large covariance between many 
pairs (i? m i„, i?max)- Although there are infinitely many mappings 



The first mapping introduced above offers a rather efficient sam- 
pling of the shear field. i? min can be fixed to the smallest observable 
distance i? m j n ,o for which shear correlation data are measured. This 
is given by the smallest separation for which galaxies do not blend, 
to allow for reliably measured shapes. It provides the largest range 
of angular scales accessible for a given !? max : The upper limit i? max 
can be varied between j? min and the maximum observed scale given 
by the data. 



2.6.2 Fixed ratio, r] = i? m in/'?ir 



const 



The second mapping has the advantage that a single filter function 
can satisfy the constraints (Eq.[9) for all scales. This makes it more 
convenient to combine different scales, e.g. to obtain the Fisher ma- 
trix, and might result in a universal optimal filter function. The ef- 
ficiency with respect to keeping & min constant is however reduced: 
A large ratio t], for which i? m i„ and i? max are close, samples only a 
small angular interval, resulting in a small signal-to-noise. A small 
j] on the other hand means that we can not go to very small scales 
with !? max : Because of the minimum observable galaxy separation 
# mi „.o, the smallest # max is min(# max ) = i? min ,o/?7- 

In both cases, we will use )? max as the argument of R E and 
denote it with the symbol as in SK07. 

In Appendix [A] we introduce a simple generalisation of this 
scheme to obtain an optimised function f + which fulfills the inte- 
gral constraints (Eq. [9} for all pairs (i? m ; n , i9 max ). However, the cor- 
responding signal-to-noise ratio is significantly lower than with the 
method used here, which was presented in Sect. 12.31 We therefore 
do not consider this option further. 

A remark about the analogy to the aperture-mass dispersion 
(M ap ) is appropriate here. (M ap > is obtained from the correlation 
function in a similar way as R E in Eq. {TJ, with integration range be- 
tween zero and twice the aperture radius. Its filter function depends 
on the two scales & (the integration variable) and 6 (the aperture 
radius). For different 6 one could define a different filter function. 
For convenience however, widely-used filters are functions of the 
ratio &/0, and therefore one functional form provides an E-/B-mode 
separation for all radii simultaneously. 



2.7 Optimisation 

The maximisation of a quantity Q, to be defined in the next sec- 
tion as signal-to-noise ratio and Fisher matrix figure of merit, is 
done as follows. For a given polynomial order N number of con- 
straints K > 2, we perfor m the maximum s earch of Q using the 
conjugate-gradient method iPress et alj i 19920 in the space of free 
coefficients a K , . . . , a N - l . At each step, we calculate s m (Eq.|19t and 
invert Eq. d20t to get the first K coefficients oq, . . . a K -\ ; from there 
we compute f + (Eq.ll4t and R E (Eq.[7j- 

We limit each of the coefficients a K , . . . , a N _i to the box 
[-10; 10]. In some cases the maximum-search fails, in particular 
if the polynomial order is high. The algorithm might run into a lo- 
cal maximum or hit the parameter boundary. To reduce the failure 
rate, we proceed as follows. 




*P (arcmin) *P (arcmin) 

Figure 1. Left panel: Signal-to-noise ratio for fixed i? m ; n = 0.2', as functions of the maximum scale T. The blue curve with circles and the cyan with triangles 
correspond to polynomial order N = 6 and N = 15, respectively. The red curve with squares shows the filter function from SK07. The black curve with crosses 
corresponds to the aperture-mass (M^ p > with *P being equal to the aperture diameter. Right panel: The comparison of He obtained from the optimised function 
T + with polynomial orders N = 6 and 15. 



If a maximum is found close to the parameter boundary we 
discard it and redo the maximisation with larger box size for a dif- 
ferent starting point. To initialise the maximisation for the first of 
a range of angular scales, we draw a number of random points, on 
the order 100, and start the maximisation with the point providing 
the largest Q. 

For subsequent angular scales, we use the information about 
the previous maximum to start the next optimisation. If the ratio 
1 = ^min/^max is kept constant when increasing & max = *F, we use 
the previous maximum coefficients a„ as new starting value for the 
maximisation process. This renders the search for the maximum 
point more efficient and more stable, since for small changes in *P 
the maximum will be close in a-space. 

For constant minimum scale i? m j n we devise a different strat- 
egy. It can be shown that in this case, Q is monotonously increas- 
ing with T: Let T+ be the function for scale which maximises 
Qi^Vi). For the subsequent scale x V i+ \ > define the function 
f + CF) = TfQ¥) for T < %, and T+(V) = otherwise. By 
construction, the resulting Q for scale "Ff+i is the same as for 
QCVi+i) = Q(fi). Therefore, we choose as new starting point for 
scale *P i+] the coefficients resulting from the decomposition of the 
function f + (W). This assures the subsequent maximum QCVj+i) 
found by the search algorithm to be larger or equal to the previ- 
ous one, at least in the limit of large N. In practice however, due to 
the finite order N, the orthogonal polynomials are not a good rep- 
resentation of f + for *F > where it is zero, and the resulting Q 
can actually decrease with increasing scale. 

We summarise the the combinations of Q, the angular depen- 
dencies and the choice of the starting point for subsequent scales in 
Table □ 



3 RESULTS 

In this section we define Q according to signal-to-noise and a 
Fisher-matrix figure of merit, respectively, which we maximise to 
find the corresponding optimised filter function t + . Before that, we 
comment on our choice of K and N, the number of constraints and 
polynomial order, respectively. 



Table 1. Overview of quantities kept fixed when varying the scales *P, and 
the corresponding method to determine the starting point for subsequent 
scales. 



Q 


Fixed quantity 


Starting point for subsequent scales 


S/N 


flmin = const. 


Previous maximum function, T + 


S/N 


= const. 


Previous maximum coefficients a„ 


FoM 


n = const. 


Previous maximum coefficients a„ 



3.1 Number of constraints and polynomial order 

We choose the minimum number of constraints K = 2 necessary 
for a finite-interval E- and B-mode decomposition, corresponding 
to the two integral constraints (Eq. |9). Adding the two continu- 
ity constraints (Eq.|16l> resulted in significantly lower values of Q. 
Some of the resulting functions T ± showed strong variations and 
narrow peaks for \x\ near unity. This indicates that the continuity 
constraint is not very "natural" but represents a strong restriction 
on the optimised filter functions. The price that has to be paid for 
continuity is then a function which fluctuates strongly, which may 
be problematic when applied to noisy data. 

Larger values of N improved Q to some extend but at the same 
time increased the occurrence of local maxima found by the search 
algorithm, reducing stability and reproducibility of the results. The 
functions f+ sometimes showed a high number of oscillations. A 
good choice for the polynomial order A' was found to be 6, equiv- 
alent to 4 free parameters a„ for K = 2. A lower A' resulted in 
significantly smaller Q. 

For the remainder of this paper we therefore choose N = 6, 
K = 2, if not indicated otherwise. 



3.2 Shear covariance and cosmology 

For the optimisation process we rely on a model shear correlation 
function and covariance. The fiducial cosmology for our model is a 
flat ACDM Universe with £l m = 0.25, Q b = 0.044, h = 0.7 and 
erg = 0.8. We use the non-linear fitting formula of ISmith et al.l 



K=2 ^min = 0.2' 



N=KK=2 fl min = 0.2' 



K=2 ti = 1/10 



K=2 x\ = 1/50 




d2003h together with the lEisenstein & Hvj {l998) transfer function 
for the matter power spectrum. The redshift di stribution of source 
galaxies is the best-fit model of iFu et al.l d2008h which has a mean 
redshift of 0.95. 

We use the covariance m atrix C++ of the shear correlation 
function £ + from lFu et al.l d2008h . corresponding to the third data re- 
lease of CFHTLS-Wide. This includes ellipticity noise and cosmic 
variance as well as the residual B-mode added in quadrature. The 
Gaus sian part of the covariance w as calculated using the method 
from lKilbinger & SchneideJ d2004h . Non-Gaussian corrections on 
small scale were applied according to Sem boloni et al.l d2007h . 

The 7?E-covariance matrix, <7?|>, is an integral over C++, 



(ftlOFi,^)) = dV&T+(ft) dtf' &'T+{&')C++{&,&'). 

(34) 

We use the two upper scale limits and 4*2 as arguments of (7? E ); 
it also depends on the two lower scale limits J? m j n> i and J? m i n ,2- 



3.3 Signal-to-noise ratio 

The first criterion for which we optimise the filter function is the 
signal-to-noise ratio 



S/NOF) 



1/2 ' 



(35) 



The variance ('Kg( v P, Y)) is the diagonal of Eq. J34t . As mentioned 
before, the signal-to-noise ratio does not depend on the normalisa- 
tion of f+ (Eq.|2TJl- 



3.3.1 Signal-to-noise for fixed i? m ; n 

The signal-to-noise is calculated as a function of T = i? raax , keeping 
& mm constant. We choose # min = 0.2' which is a typical (albeit 
conservative) lower limit where galaxies from ground-based data 



can be well separated. For each scale 4* we obtain an optimised 
filter function T+, as discussed in Sect. l2.7l (see Table[T](. 

As can be seen in the left panel of Fig. [T] the optimal filter 
function with N = 6 and K = 2 results in a much higher signal-to- 
noise than for original ring statistic function (7?/?)e, using the filter 
function Z+ from SK07. The new filter is superior to the aperture- 
mass dispersion. Note that we plot S/N for (Mj p ) as function of the 
aperture diameter instead of the radius, to have the same maximum 
shear correlation scale 4* as for He- 

Although the optimal signal-to-noise is expected to be mono- 
tonic as function of f, this is clearly not the case. However, when 
increasing the polynomial order /V to 15, the signal-to-noise is 
nearly constant for ¥ > 19' and larger than for N = 6 (see Fig.[T]l. 
As a drawback, the S/N-curve for /V = 15 is less smooth due to 
the difficulty of finding the global maximum. We did in general not 
find significantly larger values for S/N with N larger than 15. 

The similar shape of S/N for the different cases 7? E , (M^ p ) and 
(7?'??}e is an imprint of the covariance structure of the shear corre- 
lation function. The shape is modified stronger for high N, where 
the peak at around 20 arc minutes disappears. 

The shear function He is plotted in the right panel of Fig.Q] It 
has a similar shape as CRKJe from SK07, and also as (Af^ p ). This 
is reflecting the fact that all functions correspond to narrow filters 
and are band-pass convolution of the power spectrum. 

The normalised optimal filter functions for two angular scales 
are shown in the left two panels of Fig. [2] The similarity of the 
functions for different scales shows the relatively weak dependence 
of the filter shape on angular scale. 



3.3.2 Signal-to-noise for fixed t] 

Instead of fixing & mm , we now leave 77 constant and change i? min 
along with The signal-to-noise ratio increases with decreasing 77 
(left panel of Fig.[3}. This is not surprising since a larger 77 means a 
smaller range of angular scales. For 77 = 1/50 the optimal S/N ex- 
ceeds the one using the aperture-mass dispersion. The optimal filter 
functions have similar shape to the previous case of a fixed & m i n (see 
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Figure 3. Left panel: Signal-to-noise ratio as function of Y, for fixed values of rj with N = 6 and K = 2. The blue curves correspond to the optimised function 
f + . The green curve is obtained by applying the the optimal function 7 + of *P = 19' to all scales, instead of optimising each time separately, as done for the 
blue lines. The red curves shows the filter function Z + from SK07. The black curve with crosses corresponds to the aperture-mass with aperture diameter ?? max . 
Filled and open symbols represent r\ = 1/50 and 1/10, respectively. Right panel: 7?e for fixed values of r\ of 1 / 1 (open circles) and 1/50 (filled) . For the green 
line (with filled triangles), the optimal function f + for *P = 19' was used on all scales, instead of optimising each time separately as for the blue lines. 
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Figure 4. Figure of merit (FoM, Eq. I37I for the optimised filter function 
T + with rj = 1/10 and rj = 1/50 (blue curves). The red curves shows the 
FoM for the filter function Z + from SK07. The black curve with crosses 
corresponds to the aperture-mass with aperture diameter i* max . Filled and 
open symbols represent = 1/50 and 1/10, respectively. 




Figure 5. The normalised functions f + and T-, optimised for the Fisher 
matrix figure of merit (Eq. l37t with 'Pmax = 222'. Two cases, n = 1/10 and 
1/50 are shown. 



large and this case of fixed rj shows a S/N which is mostly larger 
than the aperture-mass dispersion. 

The shear function 7?e is plotted in the right panel of Fig. [3] it 
is very similar in shape as in the case of fixed i? min (Sect. l33?fl , Ta- 
ble [2] shows the polynomial coefficients of the corresponding filter 
function T+, 



the right two panels of Fig. [2]>. We use the previous maximum co- 
efficients a„ as starting point for subsequent scales as discussed in 
Sect. l2.7l (see TableQJ. Unlike in the previous case of fixed we 
do not expect S/N to be monotonous as function of *P, because i? min 
increases with V. 

Since 77 is constant, each filter function provides a valid E- 
/B-mode decomposition for any given scale. We use the filter op- 
timised for *F = 19', where the highest signal-to-noise occurs, and 
apply it to the other scales (see the green curve with triangles in 
the left panel of Fig.fJJ. As expected, the signal-to-noise for scales 
*F j= 19' is lower than in the previous case, where the optimisation 
was done for each scale individually. The difference however is not 



3.4 Fisher matrix 

To optimise the filter function we take now the alternative approach 
of minimising the errors on cosmological parameters from our new 
second-order shear statistic. To that end, we use the Fisher matrix 
F, given by 



-1 8K E C¥d dRzOVj) 
dp a dpp 



(36) 



The cosmological parameters are comprised in the vector p. As in 
the case of the signal-to-noise ratio, the Fisher matrix is indepen- 
dent of the normalisation of the filter function T + . 

The quantity Q to be maximised is the inverse area of the error 
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Figure 6. Covariance matrix (left panels) and correlation matrix (right) of Ke, optimised for the figure of merit (Eq. [37} with q = 1/50 (top), and of the 
aperture-mass dispersion with T, being equal to the aperture diameter (bottom). The colors correspond to the same levels for Re and (AfjL). In the right panels, 
the contour lines start from the innermost value of 0.9 and are spaced by 0.1. 



ellipsoid in parameter space, given b y the Fisher matrix . In two 
dimensions this figure of merit (FoM. lAlbrecht et"al1l2006h is 

FoM- 1 = jr(cniOfe - o 2 ^' 1 ; a% = (f- 1 ).. . (37) 

lEifler. Schneider & Krausd d2009h used the quadrupole moment 
determinant q of the likelihood function to quantify the size of 
the parameter confidence region. In case of a Gaussian likelihood 
(which correspond to our Fisher matrix approximation) in two di- 
mensions, the relation FoM~' = nq holds. 

We keep q constant, allowing for a single filter function T + 
to provide the E- and B-mode decomposition for each scale *P in 
Eq. d36t , and also for the covariance matrix. This requirement is 
not a necessity since the covariance between scales can be easily 



generalised to different filter functions. However, we choose this 
approach for simplicity. The cosmological parameters we consider 
are £2 m and cr 8 . 

In Fig.|4]we compare the figure of merit for the optimised filter 
function and the aperture-mass dispersion. For a given maximum 
scale "Fmax we vary in Eq. i36\ between 4.0' and f^. With q = 
1/50, the minimum angular scale is min(i? m j n ) = 4.8 arc second. 
For q = 1/10 we use the same angular scales, starting with 4.0', 
for consistency with the case q = 1/50. Alternatively, by using the 
same minimum angular scale of 4.8 arc second, the smallest ¥ can 
be in principle as small as 4.8"/^ = 0.8'. This addition of small 
scales results in a higher FoM which is comparable to the one for 
<<>• 
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Figure 7. The average correlation r (Eq. l39t for scales which are separated 
by the ratio x, for K E (solid, blue lines) and (M^ p ) (dashed, black curve). 
The error bars indicate the scatter when averaging over different scales T. 



We choose the same range of scales for the aperture-mass dis- 
persion, i.e. we vary the aperture diameter between 4.0' and '¥ mai . 
Note that although the minimum angular scale is theoretically zero, 
in practise we are limited by the smallest scale for which the £+- 
covariance matrix is calculated which is 3" in our case. 

The normalised optimal filter functions f + are shown in Fig.[5] 
They have a similar shape as the functions optimised for S/N (see 
Fig. [2}. Table|2]shows the corresponding polynomial coefficients. 



3.5 The covariance of K E 

We calculate the covariance matrix C of K E using the optimal filter 
function from the figure-of-merit maximisation at the largest scale 
^max = 222', for a constant rj = 1/50, see Sect. 13.41 As can be seen 
in Fig.|6]the covariance is diagonally-dominated, similar to the one 
of the aperture-mass dispersion. The degree of correlation is seen 
more clearly by regarding the correlation matrix 



<*£0Fi,Y 2 )> 



(38) 



see the right panels of Fig. [6] To quantify the correlation length, we 
compute the following function, 



(39) 



which is the correlation between two scales separated by the mul- 
tiplicative factor x, averaged over all x ¥. This function is shown in 
Fig.|7] Since the lines of equal r are mainly parallel to the diagonal, 
the scatter is relatively small. The correlation of K E drops off faster 
near the diagonal than the one for (M^), and shows a slightly larger 
correlation at inter mediate distances x. The covarian ce of (KK) E 
has been studied in Eifler, Schneider & Krausd 12009) and has sig- 
nificantly smaller correlation length than the one for (M ap ). 



Table 2. Coefficients a„ of the optimised function T + for (1) S/N, *P = 
19', r, = 1/50; (2) FoM, >P raax = 222', r, = 1/10; (3) FoM, >P raax = 222', r, = 
1/50. 





S/N CP = 19') 


FoM (T 


= 222') 


n 


On (n = 1/50) 


a, fa = 1/10) 


a„0?= 1/50) 





0.1197730890 


0.009877788826 


0.1239456383 


1 


-0.3881211865 


0.1061397843 


-0.3881431858 


2 


0.5212557875 


-0.4300211814 


0.5579593467 


3 


-0.3440507036 


0.5451016406 


-0.3679282338 


4 


0.2761305382 


-0.3372272549 


0.1540941993 


5 


-0.07286690971 


0.1716983151 


0.01293361618 



3.6 Numerical limits on the B-mode K B 

We calculate the B-mode Kb from Eq. (T) by using r_ obtained 
from the optimal function 7+ (see Sect. l2.4t . Theoretically, Kb van- 
ishes but there could be a residual B-mode because the E-/B-mode 
decomposition might not be perfect. For example, there could be 
numerical issues regarding the matrix inversion of Eq. J20}. Our 
values of Kb are limited by the precision of the numerical integra- 
tion of Eq. y}. Kb goes to zero for decreasing integration step size 
and we did not find evidence for a residual B-mode. For a step size 
of A# = 5 • 10~ 4 arc seconds, we find K B /K E < 2 ■ 10~ 5 for all 
angular scales. Even if a residual B-mode should be present, it is 
straightforward to make it vanish identically. One can increase the 
polynomial order N of the decomposition by one, and determine 
the corresponding coefficient a N - l such that Kb = 0. 



3.7 Dependence on cosmology and survey parameters 

The results presented in this paper have been obtained b y using a 
specifi c covariance matrix of C++, namely the one used in lFu et all 
(2008), and by choosing a specific cosmology. Here, we briefly de- 
scribe how our results change when we modify these parameters. 

First, we illustrate the dependence on the covariance matrix. 
Instead of using the full covariance, as was done in the previous 
sections, we repeat the S/N-analysis by taking the diagonal shot- 
noise component only. This noise origins from the intrinsic galaxy 
ellipticity dispersion. As expected, the S/N and FoM increase sub- 
stantially, mainly because of the missing cross-correlation between 
angular scales in the shear correlation function. For all three cases, 
K E , (Ml v ) and {KK) E , the S/N increases monotonously with *P be- 
yond the maximum scale and does not show a peak at around 20 arc 
minutes. The relative trend between the three cases stays the same. 

We change the fiducial cosmological model by increasing <x 8 
from 0.8 to 0.9. This results in an increase of S/N of a factor be- 
tween 1.3 and 1.4, which is roughly the same for K E , (Ml p ) and 
(KK} E . Increasing the mean redshift from 0.95 to 1.19 caused the 
S/N to be higher by 1 .8 to 2, in the same way for all three cases. We 
conclude that the relative difference is not dependent on cosmology 
or the redshift distribution. 

We repeat the calculation of S/N by choosing a fixed i? rain 
which is different from our standard value of 0.2 arc minutes. The 
improvement of K E over {KK) E decreases for decreasing i^min, by 
about A(S/N) = 0.1 for Ai? min = 0.1', averaged over all scales 
T. This might be because {KK} E shows less cross-correlation be- 
tween scales which leads to a larger gain when additional scales 
are included. The gain of K E with respect to (M ap ) increases 
when lowering i?mi„, as expected, since the inclusion of more 



small scales boosts the S/N. On average, the difference is 0.05 
for each 0.1 arc minute which leads to an asymptotic value of 
[S/N(«e)]/[S/N«A^»] = 1.31. 

To check the stability of the results, we add an independent, 
uniform random variable between —p and p to each of the highest 
N - K coefficients a„ after the optimum has been found. For each 
randomisation we fix oo, . . . , uk-\ as described in Sect. l2.3l to assure 
E-/B-mode separation. The S/N and FoM are very robust against 
changes in the coefficients. For both p = 0.01 and 0.1, the changes 
in S/N and FoM are on the order p and less. 



4 SUMMARY 

We have introduced a new second-order cosmic shear function 
which has the ability to separate E- and B -modes on a finite interval 
of angular scales. This function is a generalisation of the recently 
introduced "ring statistic" (SK07). Providing the second-order E- 
/B-mode shear field correlations, general filter functions are calcu- 
lated and optimised for a specific goal. In this paper, we considered 
the signal-to-noise ratio (S/N) as function of angular scale, and a 
figure of merit (FoM) based on the Fisher matrix of the cosmolog- 
ical parameters Sl m and cr 8 as optimisation criteria. 

Our method to find the optimal filter function consists in the 
following steps: 

1. Choose the polynomial order N and number of constraints 
K > 2 and define a quantity Q to be maximised (in this work 
the signal-to-noise ratio S/N and the Fisher matrix figure-of-merit 
FoM). 

2. Draw a random starting vector of coefficients a K , . . . a N _i . 

3. For m = 0. . . K — 1 calculate s m CEq.l 191. 

4. Invert the constraints matrix equation i20l to get the first K 
coefficients a . . . a K _i . 

5. Compute the filter function T + (Eqs. [TT1[T4l 

6. Calculate the shear function He (Eq. |7]l and Q (in this work 
Eqs.[35][36j. 

7. Maximise Q. At each iteration of the maximisation process, 
repeat steps 3.-6. 

We were able to improve both S/N and FoM substantially with 
respect to the SK07 ring statistic. Moreover, we obtained better 
results than for the aperture-mass dispersion <M^ p ), even though 
the latter formally extends to zero lag and includes therefore more 
small-scale power. 

We have adapted and optimised our new second-order statis- 
tic He to a specific cosmology and survey parameters such as area 
and depth. We used a smallest scale of 0.2' which corresponds to the 
smallest separation for which galaxy images can easily be separated 
using ground-based imaging data. The cosmic shear-correlation co- 
variance corresponds to the CFH TLS-Wide third data release used 
for weak cosmological lensing dFuetalj|2003) . The coefficients 
corresponding to the optimal filter functions can be found in Ta- 
ble [2] A C-program which calculates the filter functions and the 
shear statistic 7? E is freely availably 

Our specific results can be applied to other requirements, al- 
though the results will not be optimal. Alternatively, the optimisa- 
tion method described here can easily be applied to any given sur- 
vey setting. For space-based surveys, where the galaxy-blending 
confusion limit is smaller than for ground-based observations, the 
advantage of 7? E over (M% p ) is more pronounced. 

4 http : //www2 . iap . fr/users/kilbinge/decomp_eb/ 



5 OUTLOOK 

The new cosmic shear functions He and He can be applied in var- 
ious ways in the context of detecting systematics in cosmic shear 
data and for constraining cosmological parameters. 

The optimisation for signal-to-noise is useful for detecting a 
potential B-mode in the data. In this case the E-/B-mode decom- 
position serves merely as a diagnostic of the observations. A sig- 
nificant B-mode can be a sign for residuals in the PSF correction 
or non-perfect shape measurement. It might also hint to an astro- 
physically generated B-mode signal, for example from shape-shear 
correlations or shape-shape intrinsic alignment. In both cases, the 
B-mode signal is expected to be small and it is of great importance 
to obtain a clear E-/B-mode separation without mixing of modes. 

In case of a suspected astrophysical B-mode, the separation of 
the shear field into E- and B-modes might be a decisive advantage. 
If the power spectrum contains both an E- and B-mode, P lot = P E + 
Pb, both modes mix together into their Fourier- transform of P tot , 
the shear correlation function. Thus, the E- and B-mode may not be 
uniquely reconstructed from the correlation functions. The different 
astrophysical components giving rise to Pe and P% can then only 
be separated by a E-/B-mode separating filter. 

The figure-of-merit optimisation permits 7? E to be used for ef- 
ficient constraints on cosmological parameters. The reason to use 
a filtered version of the shear cor relation function instead of the 
latter directly can be manifold (Eifleretal. 200|). A filter can be 
chosen to be a narrow pass-band filter of the power spectrum and 
is therefore able to probe its local features, unlike the broad low- 
pass band function £ + . As a consequence, the correlation length is 
much smaller and the covariance matrix close to diagonal. This has 
numerical advantages in particular in the case of many data points, 
e.g. for shear tomography. Finally, higher-order moments (skew- 
ness, kurtosis, . . .) of filtered quantities are eas ier to handle than 
higher-order statistics of the (spin-2) shear field jjarvis et alj2004l ; 
ISchneider et~ai1l20051) . 

Apart from S/N and FoM maximisation, one can think of 
other, alternative quantities with respect to which the filter func- 
tion can be optimised. For example, if a model for the B-mode 
is assumed, the signal-to-noise of Kg can be optimised to facili- 
tate the possible detection of a B-mode. Further, if cosmic shear is 
combined with other probes of cosmology; the relative gain from 
weak lensing could be maximised. This can be done for specific 
goals, for example a given dark-energy parametrisation or some al- 
ternative theory of modified gravity. However, we emphasis that the 
possibilities are restricted since the optimisation is always limited 
by the information contained in the lensing power spectrum. 

In the case of shear tomography, where the shear signal 
from different redshifts is resolved (although only partially due to 
the broad lensing efficiency kernel), one can perform a redshift- 
dependent optimisation of the filter function. This is expected to 
bring further improvements: Firstly, the projection of physical onto 
angular scales varies with redshift; using a redshift-dependent filter 
function, physical scales can be sampled optimally with redshift. 
Secondly, the power spectrum changes with varying redshift; to op- 
timise the sampling of this redshift-dependent information might 
require a redshift-varying filter. 

He is beneficial in particular on small scales, where the 
aperture-mass dispersion suffers from mode-mixing. On scales less 
than a few arc minute s there is a leakage of modes of about 10% 
dKilbinger et ail [2006). Those scales contain information about 
halo structure, substructure and baryonic physics. It is difficult to 
model those effects, the use of those small scales to constrain cos- 



mological parameters is limited. On the other hand, lensing obser- 
vations on small scales will provide important constraints on the 
physical processes involved and matter properties on small scales. 
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APPENDIX A: SIMULTANEOUS OPTIMISATION FOR 
ARBITRARY SCALES 

The optimisation scheme introduced in this paper holds for a given 
ratio of minimum and maximum scale 77 = !? m i„/)? max . In this sec- 
tion we introduce a simple generalisation of the scheme to obtain an 
optimised function T+ which fulfills the integral constraints (Eq.[9j 
for all (!? m i„, # m ax)- This comes at the expense of a poor resulting 
signal-to-noise. 

If we demand the following relation to hold 

I v = J dxx v f + (x) = 0; for v= 0,1,2,3, (Al) 

then the two integral constraints 

RI + /, =R 3 I + R 2 Ii +RI 2 + I 3 =0 (A2) 

are satisfied. However, instead of two conditions we have now four 
equations (Eq. lAlt which fix K = 4 coefficients of the decomposi- 
tion. In this case, the first four matrix elements are (cf. Eq.ll7t 

f m „ = J dxx m C„(x); m = . . . 3. (A3) 

Since there are two more integrals than in the previous case, 
the resulting function has at least two more zeros. The correspond- 
ing signal-to-noise ratio is significantly lower than in the single- 
scale case; it is even smaller than the one obtained for Z + . We there- 
fore do not consider this option further. 

This paper has been typeset from a TgX/ DTr-X file prepared by the 
author. 
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